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Abstract. We discuss decays of ultra-relativistic neutrinos over cosmological distances by solving 
the decay equation in terms of its redshift dependence. We demonstrate that there are significant 
conceptual differences compared to more simplified treatments of neutrino decay. For instance, 
the maximum distance the neutrinos have traveled is limited by the Hubble length, which means 
that the common belief that longer neutrino lifetimes can be probed by longer distances does not 
apply. As a consequence, the neutrino lifetime limit from supernova 1987A cannot be exceeded by 
high-energy astrophysical neutrinos. We discuss the implications for neutrino spectra and flavor 
ratios from gamma-ray bursts as one example of extragalactic sources, using up-to-date neutrino 
flux predictions. If the observation of SN 1987A implies that vi is stable and the other mass 
eigenstates decay with rates much smaller than their current bounds, the muon track rate can be 
substantially suppressed compared to the cascade rate in the region IceCube is most sensitive to. 
In this scenario, no gamma-ray burst neutrinos may be found using muon tracks even with the full 
scale experiment, whereas reliable information on high-energy astrophysical sources can only be 
obtained from cascade measurements. As another consequence, the recently observed two cascade 
event candidates at PeV energies will not be accompanied by corresponding muon tracks. 
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1 Introduction 

While the masses and mixings of the neutrinos are already well probed, other neutrino properties 
are less understood, such as the electromagnetic properties and the neutrino lifetime. The most 
stringent phenomenological bound on neutrino lifetime comes from the observation of neutrinos 
from supernova 1987A [1, 2], which were measured in the electron flavor. Given the uncertainty 
on the supernova neutrino flux (order 50%) and the neutrino mixing parameters, it may apply 
to the mass eigenstate ui or 1^2 ■ For the sake of simplicity, we assume that the bound applies to 
vi: Ti/mi > 10^ s/eV.'^ The (model-independent) bounds on the other mass eigenstates are less 
stringent: bounds on 1^2 lifetime are imposed by solar neutrino data, yielding T2/m2 ^ 10~^ s/eV 
for decays into invisible daughters [3-5] and T2/m2 > 10~^ s/eV for decay modes with secondary i/g 
appearance [6, 7]. Furthermore, is constrained from the analysis of atmospheric and long-baseline 
neutrino data, T^/m^ > 10^^'^ s/eV [8]. More stringent bounds can be derived when specific decay 
models are assumed (see, e.g., Refs. [9-11] for an overview). For example, solar neutrinos strongly 
limit the possibility of radiative decays [12], while for Majoron decays [13, 14], explicit bounds 
can be obtained from neutrino- less double-beta decay and supernovae [15]. Another possibility is 
the decay into un-particles [16-19]. In this study, we do not consider specific decay models, but 
focus on the phenomenology of neutrino decay, given the bounds on the lifetimes above. Especially 
neutrino telescopes [20-23] are sensitive to neutrinos with an average energy and traveled distance 

''Neutrino lifetime is usually described by = Tifi/nii, where r^^o is the rest frame lifetime of the mass eigenstate 
Vi. The origin of the mass dependence is the fact that decays scale as exp[— t/(ri,o7)] — ex;p[— (Lmi)/(i5ri,o)], *-e-, 
the rest frame lifetime Tifi is boosted by 7 = i7/mi into the observer's frame, and L (baseline source-detector) and 
E (neutrino energy) are quantities related to source and experiment. 
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many orders of magnitude larger than present neutrino experiments, and may be an interesting 
approach to probe neutrino decay. 

Neutrino oscillations are typically assumed to be averaged out over cosmological distances (see 
Ref. [24] for a discussion), and the usual flavor mixing remains to be taken into account. As far as 
decays are concerned, one distinguishes between decays into products invisible to the detector, such 
as sterile neutrinos, un-particle states, Majorons, or active neutrinos strongly degraded in energy, 
and decays into visible states, i.e., active neutrino flavors. The decays can be complete, i.e., all 
unstable mass eigenstates have decayed (see, e.g., Ref. [25]), or incomplete, i.e., the decay spectral 
signature will be visible. A complete classification of complete (visible or invisible) decay scenarios 
has been performed in Ref. [26], while incomplete invisible decay has, for instance, been studied 
for active galactic nuclei (AGNs) in Ref. [27, 28]. The description of incomplete visible decays 
is, in general, more complicated [29, 30], which is why we focus on incomplete invisible decay in 
this study (IC-40 refers to the 40-string configuration of the IceCube detector). The characteristic 
energy dependence of neutrino decay may lead to easily observable imprints in the neutrino fluxes 
in that case [27, 28]. On the other hand, it has been proposed to use the flavor composition at 
the detector, such as the ratio between muon tracks and cascades, to distinguish different decay 
scenarios [19, 25, 26, 31-33]. We will discuss both options in this study, where our main focus are 
the fluxes. 

It is often believed that decays are always complete if the sources lie far enough away, which, in 
turn, allows for the test of very long neutrino lifetimes of the order of 



As a consequence, especially objects with high redshifts may be well suited to test neutrino lifetime, 
and may potentially lead to bounds stronger than the one from SN 1987A. An example, which we 
are going to use in this work, are gamma-ray bursts (GRBs), with observed redshifts as high as 
z ~ 6 — 8; see Ref. [34] for the redshift and luminosity distribution. This leads to potentially strong 
constraints via Eq. (1.1). However, note that the relationship between distance L and redshift z 
depends for z > 0.1 (or L > 360 Mpc) on the cosmological distance measure to be used, which 
we are going to address in this study; see Refs. [35, 36] for a discussion in the context of neutrino 
oscillations, and Refs. [37, 38] in the context of pseudo-Dirac neutrinos. For the neutrino flux, 
especially redshifts z ~ 1 will dominate [39], which is a consequence of the convolution of the star 
formation rate (including some redshift evolution function) and the contribution to the total flux 
scaling as 1/d^, where di is the luminosity distance. This value is already significantly beyond the 
indicated z ~ 0.1, which means that cosmological effects have to be taken into account. 

GRBs are not only an interesting class of candidate neutrino sources because they may originate 
from high redshifts, but also because the non-observation of neutrinos from GRBs [40] has started 
to constrain the gamma-ray-neutrino connection in conventional internal shock models [41, 42], 
based on the ideas in Ref. [43]. A re-calculation of the predicted neutrino fluence from the gamma- 
ray fluence has yielded an order of magnitude lower result [44] for the same burst parameters (see 
also Ref. [45]), which is a possible explanation soon going to be tested by new data. We will use 
this new nominal prediction for the quasi-diffuse flux [44], calculated with the IC-40 bursts [46], 
in this study. Note that model-specific [44, 47] and more generic [48] astrophysical uncertainties 
may imply even lower neutrino flux predictions. On the other hand, another plausible explanation 
for why nothing has been seen yet may be that neutrinos have actually (at least partially) decayed 
between source and detector. We will investigate the impact of this hypothesis on the predicted 
fluxes. 
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One of the very recent puzzles of neutrino astronomy is the potential observation of two cascade 
events at IceCube at PeV energies [49]. In particular, assuming that these are of extragalactic 
origin, an interesting question is: Why have these been observed as cascades, and no muon tracks 
have been seen? This question may become more prominent with increasing statistics. In this 
work, we demonstrate that neutrino decay can provide an answer to that question, given the current 
constraints on neutrino lifetime; for alternative explanations, see Refs. [50, 51]. We will demonstrate 
that muon tracks may be strongly suppressed in the presence of neutrino decay compared to 
(especially electromagnetic) cascades, and that cascade measurements are much more powerful to 
find astrophysical neutrinos in this case. We will use GRBs as an example, but our conclusions can 
be applied to AGNs as well. 

This study is organized as follows: We describe the decay framework in Sec. 2, which consists of the 
general framework including cosmological distances, a simplified solution, and the proper solution 
of the redshift-dependent decay equation. We then discuss the implications for GRBs in Sec. 3, 
as one possible example of extragalactic neutrino sources. In that section, we introduce several 
benchmark GRBs and the analysis techniques, discuss the impact of decay on the neutrino fluxes, 
and introduce a phenomenologically viable model for which we show the implications for fluxes 
and flavor composition. In Sec. 4 we apply this model to a realistic stacking analysis of IceCube, 
and compare the predicted neutrino fluxes to current and future bounds. Finally, we conclude in 
Sec. 5. Note that, while we use GRBs as test case, the results from this study apply to the neutrino 
propagation from other extragalactic sources as well, including cosmogenic neutrinos. 

2 Decay framework 

Here we introduce in the first subsection the general decay framework including the aspect of cos- 
mological distances. We solve the decay equation in a simplified approach in the second subsection 
to study when a cosmological framework has to be used, and we show the proper solution in the 
third subsection. 

2.1 General framework including cosmology 

Neutrino decay of the mass eigenstate vi can be described by the usual differential equation 

= -X.N. 

dt 

Here the decay rate (inverse lifetime) Aj is given by 

\ — H^il. = 
' ~ = E 

with Tj.o the rest frame lifetime boosted by 7^ = E/rrii into the observer's frame, m = mi/rifi, 
and E the neutrino energy. Since neither the rest frame lifetime Tj^o nor the lifetime Tj = jiTifi are 
directly observable, typically k^^ = Tifl/rrii [seV""*^] is quoted as the lifetime of the neutrino mass 
eigenstate Vi. 

The time t is for ultra-relativistic neutrinos typically associated with the baseline L, which is the 
distance between source and detector. Let us assume that the neutrinos under study are produced at 
a source with a redshift z. Over cosmological distances, we therefore need to express this distance 
in terms of redshift, i.e., L{z). In cosmology, however, the relationship between distance and 
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Figure 1. Light-travel distance as a function of redshift in a ACDM cosmology (solid black line). For comparison, 
we have also included the luminosity distance (dashed curve). 

redshift depends on the cosmological distance measure, which can be luminosity distance, angular 
diameter distance, proper distance, or light-travel distance; see, e.g., Ref. [52]. For neutrino decay, 
the relevant quantity is the time traveled in Eq. (2.1), which is to be associated with the light-travel 
or look-back distance. In other words, the clock which triggers neutrino decay is directly related 
to the distance the neutrino has traveled since its production. By definition. 



with Lh = c/Hq ^ 3.89 Gpc the Hubble length [53], h{z) = H (z) /Hq, and Hq the Hubble 
constant. The Hubble parameter H(z) is defined as 



assuming a flat ACDM cosmology. We choose 0^ = 0.27 and JIa = 0.73 in the following [54], 
unless explicitly noted otherwise. 

The light-travel distance L{z) is shown in Fig. 1 compared to the luminosity distance. It is limited 
by the Hubble length Lh — 3.89 Gpc, which is the horizon beyond which ultra-relativistic particles 
cannot be seen. This has interesting implications for neutrino decay: the maximum distance 
relevant for the test of decay is limited, and therefore also the maximum time the ultra-relativistic 
neutrinos have traveled. From Fig. 1 it is clear that the notion that larger distances can probe longer 
lifetimes only holds for z < 0.1 (or L < 360 Mpc), where the cosmological distance measures are 
very similar. On the other hand, for z > 1, the light-travel distance is limited. For "typical" neutrino 
peak energies higher than about 10^ TeV for GRBs and AGNs, one can estimate that the maximal 
testable lifetimes are about lO^seV ^ from Eq. (1.1). Therefore, the lifetime limit from SN 

1987A cannot be exceeded for the electron flavor (or the corresponding mass eigenstates). Although 
neutrino telescopes may also probe Kg ^ using muon tracks (the poorly constrained lifetime of u^), 
any future galactic supernova neutrino observation, which may use different flavors and neutral 
currents, will likely provide much better bounds because of the lower neutrino energies. 

The observable neutrino flux (j)a (GeV^^ cm^^ s^^) including flavor mixing and decay (neglecting 
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possible Earth attenuation effects) is given by 



</>,„+p^ {Eo, z)= Yl (^0, z) <t>l^^,^ , (2.5) 

/3=e,/i,T 

wliere is tlie neutrino ffux witliout propagation effects and Eq is tlie observed energy. 

Note tliat we simply add the neutrino and antineutrino ffuxes here, since most observables are 
not sensitive to the neutrino polarity. The transition probability P^/? for the transition 1/^ — )• 1^/3 
including ffavor mixing and decay is given by 

Paf^{Eo,z) = Y,\Uaif\U^if ^^l^^^ = J2\Uaif\Uf,^\^Di{Eo,z), (2.6) 

i m{Eo) . 

where Ni{Eo, z) is the number of neutrinos left after traveling over a distance -L(z), N[Eq) is the ini- 
tial number of neutrinos (as a function of the observed energy), and Di(EQ, z) = Ni{EQ, z)/Ni{EQ) < 
1 is the corresponding damping factor [55]. It is the solution of Eq. (2.1), which we will address 
in the next subsections. For the mixing angles, we use the best-fit values (for normal hierachy) 
from Ref. [56], i.e., sin^ 6*12 = 0.307, sin^ 6*23 = 0.398, sin^ 6*13 = 0.0245, and the CP-violating phase 
6 = 0.89 • vr. As observables, we use the muon neutrino flux (/>^ = representing muon tracks, 

the electron neutrino flux = (pue+ue.^ representing electromagnetic cascades, and the flavor ratio 
R = (f)^/{(f>e + 4>t) [57], representing the ratio between muon tracks and cascades. The flavor com- 
position at the source is computed in a self-consistent way including the cooling of the secondary 
pions, muons, and kaons, and the helicity-dependent muon decays; see, e.g., Refs. [31, 58-60]. For 
example, for GRBs, a transition from : : (j)^ = 1 : 2 : (pion decay source) to : 1 : 
(muon-damped source) is expected, which can be predicted as a function of the source parameters. 



2.2 A simplified solution of the decay equation, and a stability paradox 

Here we solve the decay equation in the usual way used in the literature, and apply it to cosmological 
distances. Although we will demonstrate that this approach does not work, it is quite illustrative 
to understand the problems and in what redshift range simplified approaches can actually be used. 

The decay equation, Eq. (2.1), has, for Aj = const., the well-known solution 

Ni{t) = Ni(t = 0) e"^'* = Ni e"^'* , (2.7) 

where Ni is the initial number of Vi. We can easily modify that equation, using Eq. (2.3), such 
that'^ 

Afi(z) = iV,e-^>^(^), (2.8) 

where Ni{z) is the number of remaining neutrinos after traveling L{z). Note that here the redshift 
is measured relative to the origin z, whereas in cosmology z = (and quantities marked with "0") 
refer to the current epoch. 

We also need to take into account the effect of the adiabatic expansion of the Universe on the 
neutrino energy: if E{z) is the energy in the production epoch and Eq is the observed energy, then 

E{z) = {1 + z)Eq. (2.9) 

''For the sake of simplicity, we ignore neutral current cascades here, which would enter as a background. 
'^We use natural units, i.e., ft = c = 1. 
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Figure 2. Comparison of the simplified {Zi) and proper {Z2) redshift dependencies. For z < 0.1 the approaches 
are identical. The gray-shaded region illustrates the typical redshift range populated by GRBs. 



With this, the decay rate in Eq. (2.2) acquires a redshift dependence, i.e.. 

Using that in Eq. (2.8), we can therefore write the damping factor needed for Eq. (2.6) as 

Kj L iz) 



A(Eo,.)=exp^^--j^J (2.11) 

with the definition of L{z) from Eq. (2.3). In order to isolate the redshift-dependent part of the 
evolution, we define the dimensionless function 

{z) = exp ^^'\ ^ , (2.12) 



Lh-{1 + z) 
such that Eq. (2.11) can be re- written as 

A(^o,2) = [^i(^)]"^ . (2.13) 

Here the exponent depends on Ki and observed energy Eq only, whereas the base depends on 
redshift only. One can easily see that for Eq — )• 00, the damping factor Di{EQ,z) — t- 1, which 
corresponds to stable neutrinos. For finite £^0 ^-nd large decay rates ^ Eq/Lh, however, one 
finds Di{EQ, z) = 1 (stable neutrinos) for Zi = 1, Di{EQ, z) — )• (complete decays) for Zi > 1, and 
Di(EQ,z) — )• 00 (exploding solution) for Zi < 1. Because large enough decay rates correspond to 
complete decays, only the second case can be a physical solution. 

Since Zi is a function of redshift, we show the functional dependence Zi{z) in Fig. 2. Obviously, 
Zi > 1, which means that the number of particles does not explode (third case excluded). However, 
for large z, Z\ — )• 1, which means that the neutrinos from high redshifts are stable, and the notion 
of complete decays does not exist. From Eq. (2.13), they are even stable independently of the 
value of Eq^ which is counter-intuitive, since more energetic neutrinos should live longer than less 
energetic ones. So what is wrong with this approach? Obviously, introducing the energy dependence 
Eq. (2.9) after the solution of Eq. (2.1) does not take into account that the neutrinos lose energy 
continuously, which changes the solution of the differential equation. That is, the energy assigned 
to the neutrinos is chosen to be the one at the production point, which is on average too high, and 
the neutrinos are hence too long-lived. While one would naively expect that this effect is small, we 
have encountered a severe contradiction here, which we are going to resolve in the next subsection. 
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2.3 Solution of the redshift-dependent decay equation 



The proper solution is, of course, to re-write Eq. (2.1) in terms of redshift using the redshift- 
dependent decay rate in Eq. (2.10): 

dNi{EQ,z) ^ KidL NiiEo,z) 

dz Eq dz 1 + z ' ^ ' ' 

Note that in the following, we choose z = as the current epoch, i.e., we integrate from z to 0. 
Therefore, from the definition of light-travel distance Eq. (2.3), we have 



dL Lh 



dz {l + z)h{z) 

The solution is 



(2.15) 



kLh dz' 
~W Jo {l + z' fh{z') 



D.,{Eo,z) = exp ( I ,^ ) = [Z2 (z)] , (2.16) 



where we have isolated the z-dependent part through"^ 

dz' 

Jo (1 + -Z ) h{z') 
Z2{z) = e'^^'\ (2.18) 

The function Z2 {z) is also shown in Fig. 2. It is now monotonously increasing with redshift, and 
larger than unity. This means that for large decay rates Hi 3> Eq/Lh, one finds Di{Eo,z) — )• 
(complete decays) for z > 0, as expected in a self-consistent framework. In addition, since the 
maximal distance is limited by the Hubble distance, the base Z2 in Eq. (2.16) is asymptotically 
limited. This has interesting implications for the notion of complete decays, which one can read 
off from Eq. (2.16): The condition for complete decays is kLh/Eq S> 1, or, equivalently. 

First of all, this confirms our earlier estimates that energies considerably lower than 1 TeV are 
needed to test lifetimes comparable to the SN 1987A bounds. Second, it is a condition on energy, 
which is independent of distance (or redshift). Therefore, over cosmological distances (for z > 1), 
"complete decays" is an energy-dependent concept, and has nothing to do with distance. 

A useful analytical approximation for Z2 is 

Z2 (z) ~ a + be-^' , (2.20) 

with^ 

C a ~ 1.71 

= -0.71 . (2.21) 




"^Note that the function l2(z) was aheady encountered in the context of oscillations, and corresponds to Eq. (6) 
in Ref. [35]. 

■^Actually, a and h are not independent, since Z2{z — 0) = 1 must be satisfied. This implies that b — 1 — a. 
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Table 1. Properties of the four bursts discussed in the following, see Ref. [61] for SB ("Standard Burst", similar 
to Refs. [43, 62[), Refs. [63, 64[ for GRB080916c, Refs. [63, 65[ for GRB090902b and Refs. [63, 66[ for GRB091024. 
The luminosity is calculated with = 47rd|, • J^^/Tgo with T-i the fluence in the energy range IkeV — 10 MeV. 
Adopted from Ref. [67]. 

Hence, the asymptotic value x's, a for z — )• 00. It depends on cosmology, as it is illustrated in 

Fig. 2 for different values of and f^A- Clearly, the higher the cosmological constant contribution 
(for a flat universe), the more efficient neutrino decays will be for large redshifts. In principle, if the 
neutrino lifetime was known from different sources, one might even use this effect to independently 
probe cosmology by comparing low and high redshift populations of certain objects - although the 
expected precision is probably not very high. 

3 Application to GRBs 

We introduce our benchmark GRBs and method in the first subsection. Then we show the impact 
of decay on the benchmark fluxes and compare the simplified and proper solutions. In the last 
subsection, we define a phenomenologically viable model and show the implications for fluxes and 
flavor ratios. 

3.1 Benchmark GRBs and method 

In the following, we apply the above decay framework to GRBs, some of which have been identified 
with high redshifts. Note, however, that our results can be equally applied to AGNs or other 
objects. We choose four benchmark bursts, for which the parameters are given in Table 1. In 
these cases, the observed photon spectrum is assumed to be a broken power law with spectral 
indices a and /3 below and above the break at e^^break) respectively. The photon fluence is J^y in 
the energy range IkeV — 10 MeV. The estimated boost factor is T, the time variability is t^, and 
the burst duration (where 90% of the counts are observed) is Tqq. As far as the selection of bursts 
is concerned, the standard burst "SB" is a burst for which most parameters, such as a, /3, z, and 
LJ^°, are frequently used in the literature. In addition, three different recent Fermi-measured GRBs 
have been chosen as examples: GRB 080916C, GRB 090902B, and GRB 091024. GRB 080916C 
has been selected, because it is one of the brightest bursts ever seen, although at a large redshift, 
and one of the best studied Fermi-LAT bursts. The gamma-ray spectrum of GRB 090902B has a 
relatively steep cutoff. GRB 091024 can be regarded as a typical example representative for many 
Fermi-GBM bursts [63], except for the long duration. Note that the middle two bursts have an 
exceptionally large F > 1000, whereas F ~ 200 for the last burst. All three observed bursts have 
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in common that the required parameters for the neutrino flux computation can be taken from the 
hterature; see the figure caption for the references. 

Once the parameters in Table 1 are given, the neutrino fluxes can be predicted, using methods such 
as in Refs. [41, 42, 44]. We use the most updated numerical method from Ref. [44] (see Ref. [39] 
for details), which has corrected earlier flux predictions down by about one order of magnitude 
by taking into account the full spectral information, the energy dependence of the mean free path 
of the protons, and the cooling of secondaries including adiabatic energy losses. In addition, t- 
channel and high-energy photo-pion production processes, neutrinos from muon, pion, kaon, and 
neutron decays, the helicity dependence of the muon decays, and flavor mixing are included. As 
one of the additional assumptions of the method, we assume that the protons carry ten times more 
energy than the electrons, and that there is energy equipartition between electrons and magnetic 
field energy, i.e., we use the same parameters as in Refs. [40, 42, 46]. With these parameters it 
is possible to estimate the maximal proton energy by balancing the acceleration rate of protons 
(acceleration efficiency 0.1 used) with the faster of the synchrotron loss and adiabatic loss rates, 
which are assumed to be the dominant energy loss rates. 

3.2 Impact of decay on neutrino fluxes, and comparison of methods 

In order to show the impact of neutrino decay on the neutrino fluxes, let us first of all choose a 
(hypothetical) model where all mass eigenstates have the same lifetimes. We show in Fig. 3 the 
predicted muon neutrino fluxf as a function of neutrino energy for the four benchmark bursts in 
Table 1 for this model. The different lifetimes used for the different curves are shown in the plot 
legend. Thick curves correspond to the proper solution in Sec. 2.3, and thin curves to the simplified 
approach in Sec. 2.2 for comparison. Note that this model is hypothetical in the sense that the 
lifetimes are chosen such that the SN 1987A bound for one of the mass eigenstates must be violated. 
However, it is the simplest model to illustrate the implications of neutrino decay. 

The actual shapes of the different predictions in the stable case ("no decay") depend on the burst 
parameters. For large enough magnetic fields, the characteristic wiggles from the secondary (pion, 
muon, kaon) cooling can be seen. For instance, for "SB" (upper left panel), the three peaks 
correspond, from left to right, to muon decays, pion decays, and kaon decays, respectively. In some 
cases, the peaks are less pronounced, especially if adiabatic cooling on protons and secondaries is 
important (such as for GRB 080916C); c/.. Fig. 13 in Ref. [60]. 

As far as the impact of neutrino lifetime is concerned, note that all bursts in Table 1 have redshifts 
larger than one. As a consequence, we can use Eq. (2.19) to estimate where the decays will be 
complete. In fact, the energy Eq, shown in the plot, below which the flux is practically switched off, 
can be computed with this equation. For instance, in the upper left panel, the flux is suppressed 
below 4000TeV for = 10^ seV~^. This condition is the more accurate, the higher the redshift 
of the burst is. For instance, for GRB 091024 (z ~ 1-1), the suppression is somewhat weaker 
(shifted to smaller energies) because the base Z2 in Fig. 2 has not yet reached its saturation value. 

Let us now compare the simplified approach in Sec. 2.2 (thin curves) with the proper solution in 
Sec. 2.3 (thick curves): Obviously the suppression in the proper treatment is stronger, which is 
expected from Z2 ^ Zi for large redshifts. The difference between the curves also increases with 
redshift, and can be quite significant for large redshifts (c/., GRB 080916C). Now one may argue 
that for the individual burst, one would measure a too short lifetime in the simplified approach, 
which just needs to be corrected. This argument, however, breaks down if different bursts at 

'Here the neutrino flux is just the neutrino fluence divided by Tgo. 
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Figure 3. Predicted muon neutrino flux as a function of neutrino energy for the four benchmark bursts in 
Table 1 (panels). The different curves in each panel correspond to no decay (stable neutrinos) and three different 
lifetimes, as given in the legend. Thick curves correspond to the proper solution in Sec. 2.3, and thin curves to the 
simplified approach in Sec. 2.2 for comparison. Here we (hypothetically) assume that all mass eigenstates have the 
same lifetime, as indicated in the legend. 



different redsliifts are compared. Let us illustrate this with one example. For the diffuse neutrino 
ffux, it is expected that the peak contribution comes from a redshift z ~ 1 [39] (Fig. 3), which means 
that low redshift (z < 1) and high redshift (z > 1) bursts contribute almost equally. In a diffuse (or 
a stacked) neutrino flux, the fluxes from these individual bursts are summed over. For z <C 1, the 
simplified and proper approaches are identical, whereas for z > 1, the properly calculated fluxes 
are more strongly suppressed for a certain lifetime. Hence, if the simplified approach is used, there 
will be a contradiction implied. 

It is, in fact, also interesting to study if z ~ 1 as the dominant contribution actually holds in 
the presence of decay. We have tested this for different values of Eq. For instance, we choose 
Eq ~ 10^ GeV, which is the peak for the SB spectrum, where the neutrino telescopes typically have 
good sensitivity. Eq. (2.19) can then be used to estimate that the flux at this energy is suppressed 
for < 400 seV-^ Fr om the convolution of the GRB rate (following the star formation rate, 
possibly with some correction), the l/d'j^ drop of E'^cf), and the decay damping, one can show that 
at this energy, the main contribution does not come from z ~ 1 anymore, but from very low z 
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(peaking at z = 0) if < 400 seV^ . This is expected, since the neutrinos from the low redshift 
bursts will not have decayed yet. If, on the other hand, > 400seV~^, the bursts from z ~ 1 
will dominate, as for stable neutrinos. 

Finally, comparing the lifetimes used for Fig. 3 with the SN 1987A bound, one can easily see that no 
strong decay suppression can be obtained if this bound is applied to all mass eigenstates. Therefore, 
the neutrino lifetimes of the different mass eigenstates have to be different for a phenomenologically 
allowed scenario, where the most plausible case, satisfying the SN 1987A constraint, might be that 
ui is stable. In fact, this is for the discussed GRB fluxes equivalent to saturating the SN 1987A 
bound on z^i. Therefore, we only discuss this case in the following, and we only use the proper 
solution of the decay equation. 

3.3 Fluxes and flavor ratios for i^i stable 

As discussed above, we focus in the following only on the case of vi stable, which satisfies the SN 
1987A constraint. Recall that the reason for this constraint is that only the electron flavor has 
been measured by the SN 1987A observation, which consists mainly of ui. On the other hand, 
we assume that f2 and 1^3 decay with the same lifetime, for the sake of simplicity, with lifetimes 
satisfying the 1^2 (and, consequently, ^3) bound. The mass eigenstate 1^2 is about an equal mixture 
among z^e, i^fi, and Ur, whereas 1^3 is a mixture between and v-j- with a small contribution of fg- 
The decay of U2 will affect the electron flavor somewhat, which however can be absorbed in the 
uncertainties of current supernova flux predictions. The decay of 1^3 will be hardly observable in 
the electron flavor. 

We show the electron, muon, and tau neutrino fluxes for this model in Fig. 4 for the standard 
burst SB and different values of the lifetime (panels). In the stable case (upper left panel) at the 
peak, where the flavor composition z^e : f/^ : i^t is approximately 1 : 2 : at the source, flavor 
equilibration at the detector is roughly achieved. Even for larger energies, the fluxes of all flavors 
are roughly comparable on a logarithmic scale - we will discuss the flavor ratio below. In the most 
extreme decay case (lower right panel), both 1^2 and 2^3 decay with a lifetime corresponding to the 
current 1^2 bound. In this case, the electron neutrino flux survives (because vi is stable), whereas 
the muon and tau neutrino fluxes are suppressed. The asymmetry between the muon and tau 
neutrino fluxes mainly comes from the large value of ^13, the chosen 6cp, and the non-maximal 
^23- As a consequence: 

• Muon tracks are suppressed by about one order of magnitude compared to the stable case. 

• Electromagnetic cascades are observable at about the same level as in the stable case. 

• Hadronic cascades are also strongly suppressed. The contributions of the leptonic tau decay 
channel into leptons can be neglected as well, since they are doubly suppressed by tau neutrino 
flux and the tau branching ratios. 

• Neutral current cascades are reduced to about 1/3 of the stable case. 

• Glashow resonances should occur at roughly the expected level as in the stable case. Note, 
however, that the right energy has to be hit (6.3 PeV), and that Glashow resonant events 
are expected to be suppressed for neutrino production from photohadronic interactions, see, 
e.g., Ref. [59] for a detailed (quantitative) discussion. 

It is therefore clear that this scenario can easily explain why only cascades have been recently 
observed at PeV energies, and will be dominating with increasing statistics, in spite of the eflective 
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Figure 4. Predicted muon neutrino flux as a function of neutrino energy for the benchmark burst SB in Table 1. 
The different panels correspond to different lifetimes, as given in the panels. The different curves correspond to 
different flavors (after flavor mixing), as given in the legend. Thin dashed curves shown the stable case (upper left 
panel) for reference. As decay model, i^i is assumed to be stable, whereas V2 and uz decay with the indicated lifetime, 
i.e., all of the panels are consistent with the SN 1987A bound. 



areas comparable between muon tracks and cascades [49] . The PeV energies are exactly at the peak 
expected from the standard burst (SB), which means that the events may come from the diffuse 
GRB flux. 

From the lower left panel of Fig. 4, one can read off that the qualitative conclusions do not change 
even if the V2 and lifetimes are increased by four orders of magnitude. In this case, only the 
rightmost peak (neutrino production from kaon decays) is not affected by the decays. A very 
interesting case are intermediate lifetimes, such as ~ lOOseV"^ in the upper right panel. Here 
the muon (and tau) neutrino flux dominates above about 1 PeV, whereas below 1 PeV the electron 
neutrino flux dominates, i.e., there is a spectral swap induced by neutrino decay. 

The implications for the flavor composition can be better seen on a linear scale. We therefore show in 
Fig. 5 the flavor ratio R (muon tracks to cascades) as a function of neutrino energy for two different 
bursts (panels) and different lifetimes (curves). One can easily see the characteristic transition from 
a pion beam source (initial flavor composition u^.v^Wr of 1 : 2 : 0, yielding R ~ 0.52) to a muon 



- 12 - 




Figure 5. Predicted flavor ratio R (muon tracks to cascades) as a function of neutrino energy for the first two 
benchmark bursts in Table 1 (panels). The different curves in each panel correspond to no decay (stable neutrinos) 
and three different lifetimes, as given in the legend. As decay model, vi is assumed to be stable, whereas 1/2 and ^3 
decay with the indicated lifetime, i.e., all of the lifetimes are consistent with the SN 1987A bound. 



damped source (initial flavor composition of : 1 : 0, yielding R ~ 0.67) at higher energies, which 
comes from the energy losses of the muons in the magnetic field. It depends, apart from particle 
physics parameters, on the value of B' , the magnitude of the magnetic field in the shock rest frame; 
see Ref. [60] for a review. The different panels show two different cases for this transition. 

It is easy to show that if only one mass eigenstate is stable, an asymptotic value of the ratio R 
is reached for complete decays which does not depend on the initial flavor composition (see, e.g., 
Refs. [26, 68]). This asymptotic value R ~ 0.14 is marked in Fig. 5, and it is always reached 
for low enough energies in the presence of decay. The low value simply reflects that muon tracks 
are strongly suppressed compared to the cascades. The transition towards this asymptotic curve 
depends on the lifetime, of course. Depending on the chosen lifetime and magnetic field in the 
source, the pion beam flavor ratio may be actually observable in a certain energy range, such as 
for = lO^seV"^ in the left panel, or not, such as for = lO^seV"^ in the left panel. 
The presence of different flavor ratios may be useful to distinguish different decay scenarios; see 
Refs. [26, 33]. 

4 Consequences for realistic GRB stacking analyses 

We now apply the scenario in Sec. 3.3 to a realistic sample of GRBs, as it has been actually used in 
a stacking analysis. Remember that from the SN 1987A constraint, we have been forced to choose 
i/i stable for all practical purposes, whereas we have assumed that 1^2 and ^3 decay with equal rates. 

We show in Fig. 6 the prediction for the quasi-diffuse muon (left panel) and electron (right panel) 
neutrino fluxes from the 117 bursts in the IC-40 stacking sample as curves labeled "no decay". The 
left curve corresponds to the nominal prediction in Fig. 3 of Ref. [44], updated with the current 
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Figure 6. Prediction of the quasi-diffuse muon (left panel) and electron (right panel) neutrino flux from the 
bursts used in the IceCube stacking analysis. The different curves in each panel correspond to no decay (stable 
neutrinos) and three different lifetimes, as given in the legend, which are applied to each burst in the stacking 
sample individually. As decay model, v\ is assumed to be stable, whereas V2 and decay with the indicated 
lifetime, i.e., all of the lifetimes are consistent with their current bounds. The limit from IC-40+59 is shown [40], 
as well as an extrapolated one for IC-86 (ten years, v4j,§*^ ~ 3 x see, e.g., Ref. [69]). For the prediction, the 

117 bursts in the IC-40 sample have been used ]46], the "no decay" corresponds to the nominal prediction in Fig. 3 
of Ref. ]44]. 



values of tlie mixing angles. It uses exactly the same bursts and parameters as the original IceCube 
analysis, but the numerical (updated) method for the flux prediction from Ref. [44]. Note that the 
quasi-diffuse flux normalization hardly depends on the burst sample used, but only the shape (c/., 
dashed curves in Fig. 1 of Ref. [40]). This means that the shape will be affected by the actual bursts 
in the sample, whereas the normalization will remain approximately constant - modulo statistical 
fluctuations from the quasi-diffuse flux extrapolation; see Refs. [39, 70]. We also show in Fig. 6, 
left panel, the current stacking limit from muon tracks, as well as an extrapolation for the full-scale 
experiment. While the current limit does not yet exceed our nominal prediction, the full scale 
experiment will finally exceed it for the chosen standard values of the astrophysical parameters. 
Note that in the original IceCube analysis, z = 2 has been assigned to a burst if the redshift has 
not been measured, while the pion production efficiency has been computed (for long bursts) with 
Ll^° = lO^^ergs"^. Since we stick to these rules, most bursts will have a redshift z > 1, and 
cosmological effects will be important for decay. 

In Fig. 6, we also show the flux predictions for different neutrino lifetimes. In the extreme cases, 
< lseV~^, the muon neutrino predicted flux is suppressed to below the expected limit from 
IC-86 over ten years. However, in these substantial electron neutrino flux is expected; see 

right panel. It can be read off from the figure, that gamma-ray burst neutrinos should be detected 
with the full scale experiment in cascades if the effective area for cascades is at least 25% of the one 
for muon tracks at 1 PeV.^ On the other hand, no reliable information on astrophysical neutrino 

^Unfortunately, there are no GRB cascade analyses available up to this point, which means that reliable limits 
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sources, such as GRBs, can be obtained from muon tracks only. In particular, it is not clear if no 
neutrinos are found because they decay, or because the baryonic loading in GRBs is smaller than 
anticipated. 

A very interesting case is the intermediate value k^^ ~ 100 seV^^ and its implications for muon 
tracks; see left panel. In this case, only the low energy component of the flux is damped, and 
only neutrinos above a certain threshold energy are observable. It is expected that the full scale 
experiment may detect neutrinos in such a case, but the interplay with the (energy-dependent) 
detector response becomes important (see Ref. [71] for a more detailed discussion). For instance, 
better information may be obtained from down-going than up-going muon tracks, because these 
exhibit a better sensitivity at higher energies (c/.. Fig. 3 in Ref. [71]). As a consequence, the 
sensitivity may be dominated by a different viewing window and therefore different bursts compared 
to the "no decay" case, which may have different characteristics. One has to keep that in mind if 
one compares the sensitivities for these cases. 

5 Summary and discussion 

We have discussed neutrino decays of ultra-relativistic neutrinos over cosmological distances, i.e., 
for redshifts z > 1. These neutrinos may originate from extragalactic cosmic accelerators, such as 
AGNs or GRBs. As class of decay models, we have considered decays into secondaries invisible 
to the detector, such as Majorons, unparticles, sterile neutrinos, or neutrinos strongly degraded in 
energy. The decays have been allowed to be incomplete, i.e., the spectral signature of the decay 
may be seen. 

We have demonstrated that simplified approaches to neutrino decay cannot be used over such 
distances, and that the differential decay equation has to be carefully solved as a function of redshift. 
We have also given a simple parameterization for the solution, which can be easily used. From a 
comparison of different approaches, we have demonstrated that conventional approaches to neutrino 
decay can be safely used for z < 0.1 (or distances L < 360 Mpc), whereas the light-travel (or look- 
back) distance triggers the clock for neutrino decays over cosmological distances. As a consequence, 
the horizon is limited by the Hubble length ~ 3.89 Gpc, which is the maximal distance scale testable 
by decay. This has several interesting consequences: First of all, the common notion "the longer 
the distance, the longer lifetimes can be tested" does not hold anymore. Second, the concept of 
complete decays (all neutrinos have decayed) cannot be assigned to long enough distances, but only 
to low enough energies, i.e., it is an energy-dependent and not distance-dependent statement. And 
third, the neutrino lifetime bound from SN 1987A, which emitted neutrinos at much lower energies, 
cannot be exceeded by high-energy astrophysical neutrinos from AGNs and GRBs tested in neutrino 
telescopes. The same argument applies to the cosmogenic neutrinos from cosmic microwave (and 
infrared) background interactions, which are at even higher energies. Note that the SN 1987A 
neutrinos have only been measured in the electron flavor, which means that this limit does not 
apply to all mass eigenstates. However, it is expected that the neutrinos from any future galactic 
supernova explosion will be measured in different flavors and neutral currents, which means that 
in long terms the best lifetime bounds are expected to come from supernova neutrinos. 

While our results can be easily applied to AGNs as well, we have used GRBs as test case. Because 
of the SN 1987A constraint, we have been forced to choose ui stable and 1/2 and ^3 decaying with 
equal rates. We have shown that, safely within the current bounds for 1^2 and 1/3, the muon 

in the right panel of Fig. 6 cannot be shown. 

''We have demonstrated that saturating the SN 1987A bound for the Vi hfetime does not change the results, 
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neutrino flux can be substantially reduced by about one order of magnitude, whereas the electron 
neutrino flux is hardly affected. This has very interesting implications for neutrino telescopes: 
First, improved bounds on the neutrino flux from muon tracks can be interpreted in terms of the 
astrophysical parameters or the possibility that neutrinos are unstable. On the other hand, the 
electromagnetic cascade bounds are hardly affected by the neutrino lifetime, which means that 
reliable conclusions on the astrophysical objects can only come from cascades, and highlights the 
importance of dedicated cascade analyses for GRBs and other classes of objects. This would be 
different, of course, if muon tracks from astrophysical neutrinos were actually observed, since their 
discovery would constrain this possibility. 

For GRBs, in particular, we have used a realistic sample of IceCube GRBs and the state-of-the-art 
technology for the prediction of the quasi-diffuse neutrino flux to illustrate the consequences. We 
have demonstrated that, for the discussed decay scenarios and our nominal flux prediction using 
the same parameters as IceCube, even the full scale IceCube experiment may not find neutrinos 
from GRBs in muon tracks after ten years, whereas the perspectives for a detection in cascades are 
actually very good. 

Finally, note that recently two cascade neutrino event candidates at PeV energies have been found. 
Our scenario provides a plausible explanation why only cascades are seen, whatever the origin of 
the neutrinos is, whereas there are no accompanying muon tracks in spite of comparable effective 
areas. For GRBs, in particular, the PeV energies are exactly where we predict the maximum of the 
quasi-diffuse flux; see Fig. 6, right panel. On the other hand, for the relevant search time window, 
we only expect 0.07 electromagnetic cascade events for the nominal ("no decay") prediction from 
GRBs. The observed two events might in this case come from a strong statistical fluctuation, or a 
significant deviation of astrophysical parameters from their assumed mean values. 
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